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Abstract. We suggest that the energy source of the observed diffuse gamma-ray emission from the direction of the Galactic 
center is the Galactic black hole Sgr A*, which becomes active when a star is captured at a rate of ~ 10~ 5 yr -1 . Subsequently 
the star is tidally disrupted and its matter is accreted into the black hole. During the active phase relativistic protons with a 
characteristic energy ~ 6 x 10 52 erg per capture are ejected. Over 90% of these relativistic protons disappear due to proton- 
proton collisions on a timescale t pp ~ 10 4 years in the small central bulge region with radius ~ 50 pc within Sgr A*, where 
the density is > 10 3 cm" 3 . The gamma-ray intensity, which results from the decay of neutral pions produced by proton-proton 
collisions, decreases according to e~' lT i>i', where t is the time after last stellar capture. Less than 5% of relativistic protons escaped 
from the central bulge region can survive and maintain their energy for > 10 7 years due to much lower gas density outside, 
where the gas density can drop to ~ 1 cm 4 . They can diffuse to a ~ 500 pc region before disappearing due to proton-proton 
collisions. The observed diffuse GeV gamma-rays resulting from the decay of neutral pions produced via collision between 
these escaped protons and the gas in this region is expected to be insensitive to time in the multi-injection model with the 
characteristic injection rate of 10~ 5 yr -1 . Our model calculated GeV and 511 keV gamma-ray intensities are consistent with 
the observed results of EGRET and INTEGRAL, however, our calculated inflight annihilation rate cannot produce sufficient 
intensity to explain the COMPTEL data. 
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£ ; 1. Introduction 



A supermassive black hole known as Sgr A*, with a mass of 
~ 3.6 x 1O 6 M (for a recent review of the properties of the 
Galactic center (GC) black hole see iGenzel & Karasl I2007I) 
is located at the GC. In addition, th ere are many high en - 
ergy sources harbored in this region dMelia & Fa lcke 2001). 
For high energy radiation ranging f rom 511 keV annihila tion 
lines detected by INTEGRAL (e.g. IChurazov et al.ll2005l) . in 
the range 1 MeV-3 MeV, flux was measured by COMPTEL 



(Strong etal. 2000 ), 30 MeV- 10 Ge V photons by EGRET 



(e.g. Maver-Hasselwander et al. 19981). an d TeV photons de- 
tected by Whipple dKosack et al.1 12004. by CANGAROO 
dTsuchiva et alj|2004 . and by HESS dAharonian et al.|[2004 
are observed from the direction of the GC. 

The spatial distribution of these four energy bands is as fol- 
lows: The 511 keV annihilation lines are emitted from a non- 
spherical symmetric extended region with about 6-8 degrees 
FWHM centered at the GC. The emission appears to be dif- 
fused and does not show any clear point source in the emission 



region. The annihilation flux is about ~ 10 43 ph s _1 or ~ 10 37 
erg s _1 . The MeV emission does not show any strong concen- 
tration in the direction of the GC and the emission is rather dif- 
fuse. An excess (about a factor of two above the prediction of 
the standard model) is found in this energy range whose origin 
is still unclear. The GeV source detected by EGRET is known 
as 3EG J 1746-2851, which has an emission region around 0.5 
degrees in radius with the flux ~ 10 37 erg s , and the GC is 
at the rim of its emission region. This region is surrounded 
by an extended region of diffuse emission whose flux is about 
1 The emission position of TeV photons with a 



10 38 erg s 
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flux ~ 10 erg s can be determined to be less than 10 pc and 
its center almost overlaps with the Galactic black hole. 

Although the emission regions of GeV and TeV photons 
have some overlap, it is unclear if they are related. In general, 
511 ke V photons are not considered to be related to both GeV 
and TeV photons. Furthermore, the annihilation line is consid- 
ered to be the most difficult problem to be explained. For the 
secondary origin of positron produced by GeV protons with 
nuclei of background gas this implies that the energy content 
in primary relativistic protons necessary to create the observed 
annihilation flux ~ 10 43 ph s _1 is about ~ 3 x 10 54 erg. This 
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amount of energy is very difficult for any known mechanisms, 
except gamma-ray bursts, occurring a t the cosmological dis- 
tance to provide. Parizot et al. ( 2005b argued that indeed the 
observed positrons result from the hypernova explosion, which 
is the progenitor of the gamma-ray bursts. However, the energy 
claim that gamma-ray bursts can reach 3x1 54 erg is assumed to 
be isotropic emission. It is generally accepted that the emission 
of gamma-ray bursts is beaming and hence the energy emit- 
ted is actually two to three order of magnitudes lower than the 

isotr opic cas e. 

IPrantzosI d2006h suggests that most positrons are produced 
in the disk but they are transported to the bulge by the regular 
magnetic field. If this is true, then similar propagation should 
be assumed for relativistic protons and electrons which also 
propagate by diffusion along magnetic field line. Their propa- 
gation in the perpendicular direction is due to random fluctu- 
ations (spaghetti-like structure). If Prantzos's model is correct 
we would observe extremely high fluxes of radio and gamma- 
ray emission from the GC. In fact the actual g alactic mag- 
netic field is quite complicated (Han et al.l 120061) . and as fol- 
lows from radio data (see Beckl 20071) thederived structure of 
the Galactic magnetic field does not c o rres pond to Prantzos's 
assumptions. iMastichiadis & Ozernovl (119941) argued that the 
gamma-rays originating from the Galactic black hole may pos- 
sibly be produced from relativistic particles accelerated by a 
shock in the accreting plasma. At the same time, the gamma- 
rays could also come from some extended feature s such as ra - 
dio arcs, whe r e rel ativistic particles are present dPohllll997l) . 
Markoff et al.l (119971) discussed in detail the gamma-ray spec- 
trum of GC produced by synchrotron, inverse Compton scat- 
tering, and mesonic decay resulting from the interaction of 
relativistic protons with hydrogen accreting onto a point-like 
sources (e.g. the massive black hole). However, the above 
models cannot produce the hard gamma-ray spectrum with a 
sharp turnover a t a few GeV, which is o bserved for the GC 
source. Recently. lOka & Manmotol (120031) have suggested that 
the gamma-rays produced in the inner portion of accretion 
disk through the decay of neutral pions created by p-p colli- 
sions may contribute to the gamma-rays observed by EGRET. 
However, their model predicted gamma-ray intensity is at least 
tw o orders of magnitud e lower than the observed intensity 



Cheng et a i1 d2006h have suggested when a star is captured 



by a supermassive black hole at the GC, the star will be ac- 
creted into the black hole and a jet may be emitted during the 
accretion process. This idea comes from the fact that accret- 
ing black holes systems are seen to be accompanied with jet 
emission. One example of such a system is the microquasars. 
Studies on microquasars reveal that these objects behave very 
differently in their high/low state. In their low state, there is ev- 
idence of jet emission, although the b ulk Lorentz fac tor of the 
jet is likely to be less than 2 (e.g. see lGallo et al.ll2003l) . When 
they are in the high/soft state, t here is evidence that the jet for- 
matio n is greatly suppressed ( Fender et al] 1 19991 iGallo et al 
20031) . However, in their "very high" state, the jet reappears 



Unlike the jet seen in the low state, the jet is very powerful and 
highly relativistic in the "very high" state of a microquasar (e.g. 
lFenderl2003l:lFender & Maccaronel2004l) . This example shows 
that it is indeed possible that the transient accreting black holes 



are accompanied by jet emission. If jet emission also occurs 
during a tidal disruption event, the ejected jet will interact 
with the interstella r medium (ISM) and decelerate accordingly. 
Wong et al.l (120071) studied the electromagnetic radiation from 
the jet produced from a star captured by the black hole. They 
compared the X-ray and optical data from some nearby galax- 
ies, which are suggested to have the recent capture events. They 
concluded that a capture event with a characteristic jet energy 
10 52 erg is capable of explaining the observed time dependent 
data. However, it is extremely difficult to observe the emission 
in radio because the synchr otron self-absorptio n has strongly 
suppressed the radio waves. ICheng et al. (2006) have assumed 
that the jet should consist of relativistic protons, which will 
gradually diffuse to a large distance away from the black hole. 
The proton-proton collisions can produce enough positrons to 
explain the observed annihilation flux of positrons from the di- 
rection of the GC. 

Below we analyze the model of central black hole presented 
in lCheng et al. ( 2006b when positrons are ejected due to stellar 
capture. Since many parameters of this process are unknown 
we shall try to estimate them from observed parameters of 
gamma-ray emission from the central region. For example in 
the framework of this model it gives a possibility to estimate 
the necessary energy release in capture processes. Our calcula- 
tion in this paper is based on a multi-capture model, in which 
the total energy release due to multi captures in 10 7 years is a 
function of the capture frequency and masses of captured stars. 

The paper is arranged as follows. In Sect. 2, we describe 
the model, in which relativistic protons will be ejected by the 
Galactic black hole when a star is captured and we introduce 
important timescales in the model. We estimate how many 
gamma-rays will be emitted from the central compact region 
and from the larger low density region. In Sect. 3, we outline 
the model calculations. In Sect. 4, we present our numerical 
results and compare our model results with the observed high 
energy radiation data from GC. In Sect. 5, we present a brief 
discussion. 

2. Model description 

The rate at which a massive black hole in a dense star cluster 
tidall y disrupts and swallows stars has been studied extensively 
(e.g. iHillsl [l975l iBahcall & Wolfl Il976l iLightman & Shapiro! 
1977b . Basically when a star trajectory happens to be suffi- 
ciently close to a massive black hole, the star would be cap- 
tured and eventually disrupted by tidal forces. After a dy- 
namical time-scale (orbital time-scale), the debris of a tidally 
disrupted star will form a transient accretion disk around the 
massive black hole, with a radius typically comparable to the 
tidal capture radius dReesll 19881) . Rees also argued that most of 
the debris material is swallowed by a black hole with a mass 
~ 1O 6 M on a timescale of ~ 1 yr for a thick hot ring, or 
~ 10 2 years for a thin cool disk. The more quantitative de- 
scription will be given later. The capture rate is essentially a 
problem of loss-cone diffusion-diffusion in angular momentum 
rather tha n energy. By a s sumin g a Salpeter mass function for 
the stars, ISyer & Ulmerl d 1999b estimated the capture rate in 
our Galaxy as ~ 4.8 x 10~ 5 yr -1 for main sequence stars and 
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8.5 x 10 6 yr 1 for red giant stars, respectively. On the other 



hand, Ma gorrian & Tremaine ill 9991) used the dynamical mod 
els of real galaxies by taking into account the refilling of the 
loss cone of stars on disruptable orbits by two-body relaxation 
and tidal forces in non-spherical galaxies, but they obtained a 
higher capture rate ~ 10 yr -1 . Therefore, the actual capture 
rate is sensitively dependent on the assumed mass function of 
stars, the stellar evolution model used, the radius and mass of 
the captured star, the black hole mass and the internal disper- 
sion velocity of stars (v ? ) around the black hole . For example, 
based on the theory of ICohn & Kulsrudl dl978l) . ICheng & Lu 
(2001) have shown that the capture rate of stars by the massive 
black hole is proportional to M^ 33 nJ- 6 v s ~ 5 - 76 , where Mbh, and 
n„ are the mass of the black hole and the star density in the 
star cluster around the black hole respectively. They obtained 
a longer capture time ~ 10 6 years, by taking v s = 10 2 km/s 
and M\,h - 3.6 x 1O 6 M . Therefore the capture time for a 
main sequence star with mass ~ 1M could range from several 
tens of thousands of years to several hundreds of thousands 
of years. It is very important to note that the correct predic- 
tion of capture rate is a ver y difficult task. Based on the ob- 



servati ons of nearby galaxi e s.lFerrarese (2002), Geb hardt et aL 
( 2002 ) and Tremaine et alj d2002h have given some simple re- 



lations between the black hole mass and the velocity disper- 



sion as Mbh 



4.58±0.52 



and Mbh 



4.02±0.32 



respectively. If 



we substitute these simple relations into the formula derived 
by Cheng & Lu, it produ c es a t otally unreasonable capture 
rate. However. [Barth etal] (2004) have shown that such sim- 
ple formulae would break down in the mass scale like the black 
hole in GC. Although the stellar capture rate is difficult to be 
determined theoretically, there are five X-ray flare events ob- 
served in nearby normal galaxies, which are believed to be 
the consequences of the stellar ca pture jDonlev et al. 2002; 
Halp ern et alj |2004; Ko mossa 2006). Based on these observed 
events the average capture rate per galaxy is about v cap ~ 10" 5 
yr 1 (below T cap = l/v cap ). 

When a star comes within the capture radius, which is given 

by 



1.4 x 10 13 Ml n m; U3 r„ cm, 



(1) 



where m, = MJM Q , M 6 = M m /1 6 Mq r, - RJR a , the sta r 



will be captured by the black hole dReesil 1 9881; IPhinnevll 19891) . 
According to the theoretical predictions, the flare results from 
the rapid release of gravitational energy as the matter from the 
disrupted star plummets t oward the black hole. For t> t pea k, the 
accretion rate evolves as, dReesil 1 9881: |PhinnevMl989h . 



1 M, 
M - 

3 tmin \ tfnin 



-5/3 



where M„ and R* are the mass and 
captured star, respectively and t pea k 

~m ) \^ ) \ io 6 m g j y r * s tne characteristic t ime for the 



the radius of 

1 .59f<m' n , t m i n 



(2) 



the 



debris to return to the pericenter dLu et alj 120061) . Recently 



Yuan et al. (2002) studied jets emission from Sgr A*. They sug- 
gested that the Chrandra observed features of X-rays from the 
vicinity of Sgr A* can be explained in terms of a coupled jet 
plus accretion disk model. The observed radiation is mainly 



emitted by the electrons in ADAF disk and in the jet. However, 
it is not known if the energy of the jet will be carried away by 
protons or electrons. Since the inertia of protons is much larger 
than that of electrons, it is logical to assume that the energy 
of the jet is mainly carried by protons. The energy distribution 
of the protons in the jet is usuall y assumed to be a power law 
and the index is taken to be 2-3 dYuanl2007l) . lYuan et ail d2005l) 
studied how much mass is carried away by the jet in the black 
hol e system, and concl uded that it is typically 1%-10%. 

Cheng et al. (2006) estimated that the maximum accretion 



energy carried away by relativistic protons is given by 



AE p ~6x KPOfo/lO-'XM./Mo) erg, 



(3) 



where rj p is the conversion efficiency from accretion power 
(Mc 2 ) into the the energy of jet motion. 

A number of timescales are important for our model; the 
proton-proton collision time scale 

t pp = {no-ppcT 1 ~ 3 x 10 V 1 yr, (4) 
the diffusion time scale 

T d = d 2 /6D ~ 10V/500 pc) 2 (D/10 27 cmV 1 ) yr, (5) 

and the ionization cooling time scale for the relativistic charged 
particle 



t coo1 = 10 8 (£/GeV)n- 1 yr. 



(6) 



Here n is the gas density, cr pp is the p-p collision cross-section, 
d is the distance from the source, D is the diffusion coefficient 
and E is the energy of charged particle. 

It is very difficult to explain the extended spatial distri- 
bution of the annihilation emission around the GC unless the 
sources of positrons are more or less uniformly distributed 
in the bulge as there is no such problem of pr opagation (e.g. 
IWang et al .1 12006; Wei denspointner et alj |2006). The problem 
is that we do not know how positrons with energies below 
100 MeV propagate through the interstellar space. From ob- 
servations, we can derive average values of the diffusion co- 
efficient which, in principle, differ from each other depend- 
ing on the analyzed spatial r e gion. O ne can find these es- 
timations in iBerezinskii et alj dl990l) ; IStrong & Moskalenkol 



(1998); Strong_etalJ(2000). Thus, from radio and gamma-ray 
emission from the Galactic halo whose semithickness is about 
several kpc, one can show that cosmic rays with energies above 
100 MeV propagate by diffusion with the coefficient in scales 
of the Galactic halo of about D h ~ 3 ■ 10 28 - 10 29 cm 2 s"'. From 
the cosmic ray chemical composition which is determined by 
particle propagation inside the gaseous disk one can find the 



value about D„ 



3 • 10 27 - 3 • 10 28 cm 2 s 



smaller than in 



the halo. Similar values for the diffusion in the local Galactic 
medium were derived from the anisotropy of high energy cos- 
mic rays emitted by nearby supernova shells. Besides, these 
estimations are strongly dependent on whether the values of D 
are spatially or energy dependent as well as if there are other 
mechanisms of cosmic ray transport in the Galaxy. 

It is unclear, of course, whether the average characteristics 
of cosmic ray propagation in the disk or in the halo can be 
extrapolated onto the region of the Galactic bulge in order to 
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describe distribution of positrons with much smaller energies. 
However, rather simple estimates for the spatial diffusion coef- 
ficie nt of Me V p ositrons give only the value, D ~ 10 27 cm 2 s 4 
(e.g. lJean et al. 2006). In the framework of our model we can 
argue that positrons ejected from a central region of the GC 
should propagate over the distance about several hundred pc 
during the time of their thermalization (~10 million years in 
order to satisfy the observations). 

The gas density di stribution in the GC is complicated. 
According to Jean et al. ( 2006h . the bulge region inside the ra- 
dius ~ 230 pc and height 45 pc contains 7 x 1O 7 M . A total 
of 90% of this mass is trapped in small high density clouds (as 
high as 10 3 cm 4 ) while the remaining 10% is homogeneously 
distributed with the average density ~ 10 cm 4 . In the 500 pc 
region the average density will drop to ~ (1 — 3) cm 4 . 

When relativistic protons are ejected after the stellar cap- 
ture, they take ~ 10 5 years to leave the central high den- 
sity region (~ 50 pc) if the diffusion coefficient of protons 
near the GC is about 10 27 cm 2 s 4 i.e. as in the Galactic Disk 



(Ber ezinskii et al. 1990). The proton collision timescale is only 
3 x 10 4 years. Therefore, most proton energy will be con- 
verted into pions, which quickly decay to photons, electrons 
and positrons. 

The most plausible observed gamma-ray intensity must be 
of the order of 



L r (50 pc) 



rj„AE p , 

—p 'cap/'pp 



6 x lO^ifc/lO^Xifo/lO-'XM./Mo) erg/s, (7) 



where T] n is the conversion efficiency from protons to neutral 
pions. 

The clumpy structure of the gas in the central 50 pc region 
is essential for estimates of the lifetime of primary protons be- 
cause the primary protons may disappear after the first penetra- 
tion into dense clouds if the gas density there is high enough. 
However, since the fraction of protons leaving the high den- 
sity region is e^ {T ^ T " ) ~ e 4 ~ 0.05 then the number of clouds 
should be large enough and the approximation of average den- 
sity is completely acceptable. 

These 5% of primary protons leaving the central region can 
propagate by diffusion through the 1 cm 4 interstellar gas to the 
distance about 500 pc. Since the p-p collision time in 500 pc is 
10 7 years, the proton injection rate into this region is almost 
constant in comparing with the diffusion timescale. Therefore 
the gamma-ray emission intensity from 500 pc region is almost 
constant and the gamma-ray power is given by 



L r (500pc) * 0.05 



r] n AE p 



* 3 x lO^Ofr/lO^XVlO^XM 'M ) erg/s. (8) 

It is very important to note that most of positrons are cre- 
ated in this high density region. However, they may propa- 
gate mainly through the intercloud medium because of the 
screening eff ect due to MHD waves excited near dense molec- 
ular clouds jSkilling & Strong! 1 19761 iDogel & Sharovl 1 19851; 
IPadoan & Scaloll2005l) . It is unclear at which precise energies 
of positrons this effect is essential but if these waves are excited 



by subrelativistic nuclei then positrons with energies as high as 
(and below) 3 00-600 MeV cannot penetrate into the clouds 
( Morfill|[l982 ). This means that a significant part of secondary 
positrons even with relativistic energies is cooled down in the 
intercloud medium only. The cooling time scale of positrons is 
as large as ~ 10 7 years. 

The expected positron annihilation rate in this case is given 

by 



Ni 



AE n 



2 x 10 43 (i7-/10 _1 )(Af »/M )s 4 , 



(9) 



The estimates (0 and © a re completely consi stent with 



the EGRET data (Maver-Hasselwander et al. 1998) and the 



estimate (|9]l w ith th e observed results of INTEGRAL (e.g. 
Churazov et alll2005h . 



3. Model calculations 

The detailed description of model calculations was given in 
Cheng et al. ( 2006i) . Here we summarize the calculation pro- 
cedure as follows. 

Since most accretion energy will be released over a very 
short timescale, we take the source function of protons as 



Q(r,E p ,t) = A(E p )6{r)6(i), 



where A{E„) oc 



(10) 



for the power-law momentum 



J P' (£j+2M„c 2 £ p )<ra+i>/2 

injection spec trum and the spectral in dex jq is taken to be be- 
tween 2-3 (cf. lBerezinskii et al.ll 19901 and reference hereafter). 

Since the acceleration processes in the jet are still unclear, 
it is not easy to determine the maximum energy of protons in 
the jet. According to the T eV observation s from AGNs (e.g. 
Mkn 421 and Mkn 501, cf. IWeekes1l2004 . if we assume that 
these TeV photons result from proton-proton collisions, the en- 
ergy of protons in the jet must be at least over 10 TeV. In fact, 
TeV photons have also been observed from the vicinity of Sgr 



A* (Kosack et al. 2004; Tsuchiva et al. 2004; Aharonian et al 



2004), it has been suggested that these TeV photons result from 
the p-p collisions and the rel ativistic protons are ejected fro m 
the Galactic b lack hole (e.g. lAharonian & Neronov 2005 alibi 
iLu et al.l2006l) . However, the exact values of the maximum pro- 
tons are not important in our problem. In order to fit the EGRET 
data, the spectral index of the proton spectrum is required to be 
close to 3. Therefore most proton energy is at ~GeV. 

The spatial distribution of the protons can easily be derived 
from the well-known equation of cosmic ray propagation, 

dn p d IdE \ n„ 

_Z_V( D V„, ) + -(-n p ) + ^ = e( r,^0. (ID 



Here D is the diffusion coefficient, 



dE 
dt 



2ne A n I m 2 c 2 W„, 
■In 



mc/3(E) \ 4ne 2 ti 2 n 



(12) 



is the rate of ionization losses and t pp is the characteristic time 



of p-p collisions. 
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Neutral pions are produced in p-p collisions and they will 
decay almost immediately to high energy photons. The emis- 
sivity of the photons produced by 7r°-decay (q Y ) can be calcu- 
lated from corresponding references for this and other equa- 
tions used for calculation s of gamma-ray an d secondary parti- 
cles which one can find in lCheng et al.1 (12006): 



s = E/m e c 2 is the dimensionless photon energy. The energy 
spectrum of in-flight annihilation 



-/ 



N a „(s) = q a „{jQ, s)n e (s)dy 



(17) 



q r (E y , t) = 2 



CO 

I 



q n {E„, t) 
yjE% - m 2 n c 



zdEyr 



(13) 



where q„ = JT N(Ep)v p nndoiE, Ep), dcr(E,E p ) is the differ- 
ential cross-section for pions, E y and E K are the energy of the 
emitted photon and the decaying pion respectively, = 

Inelastic p-p collisions produce two charged pions for 
every neutral pion. These charged pions quickly decay into 
muons, which in turn decay into positrons and electrons, with 
a resulting emissivity 



q e (E e ) = n H - 



mlc 



r dp r dE n r 

I dEp— — I I dE p n 

J dE e J B n E„ Je„,(e„) 



dcr(E„, E p ) 
dE n 



We use the generalized Fokker-Planck equation to calculate 
the positron distribution function / written in dimensionless 
variables p - pi vmkT and t = vof is 



d£ + vn(cr an + cr ce ) 
dt v 



q e (p, r) 

j_d_~ 
p 2 dp 



A(p) d -M + B( P )f 
op 



(15) 



where (cr an ) and (cr ce ) are the cross-sections for inflight an- 
nihilation and charge-exchange processes respectively, q e de- 
scribes the distribution of sources emitting fast positrons, 



vo 



^lm t .kT x 



Mp) 



\ dt h 

tdp\ 

\dth 



and B(p) 



1 [ \ at I ion \dt) synIC \at Ibrem 

The high energy gamma-rays (> 30 MeV) are mainly pro- 
duced by the decay of neutral pions, whereas the lower energy 
gamma-rays (< 30 MeV) are produced by in-flight annihila- 
tions of positrons. The in-flight differential spectrum of the 
y-rays produced by annihilation of a positron on the ambient 
electrons with density n e : 



nr 2 cn e 



7+ + 1 



2 - + ■ 
\e y + + 1 



r+ + l - e 
l \ /l 



+ 

s y+ + l 



(16) 



where y + = E + /m e c 2 is the Lorentz-factor of the positron, 
p + = sjy 2 - 1 is the dimensionless momentum of positron and 



The lower integration limit can be obtained from the kinetic 
equations and equals 



y(s) 



e 2 + (e-1) 2 
2s- 1 



(18) 



4. Numerical calculations 

4.1. Energy input constraints of multi-injection model 

From Eq. (5), we see that the energy carried away by the rela- 
tivistic protons can vary from capture to capture depending on 
the mass of the captured star and the efficiency of converting 
accretion energy into outflow of relativistic protons. It ranges 
from 10 52 erg for t] p ~ 0.01 and M*/M ~ 1 to 10 54 erg for 
T] p ~ 0.1 and M*/M Q ~ 10. In numerical calculations, we can 
only assume a mean energy injection in Eq. (10). To constrain 
the injection energy, we used the observed electron/positron 
annihilation intensity to constrain the injection energy as fol- 
lows: In p - p collisions only 4% of initial energy of proton 
is transferred to secondary electron which gives the energy of 
positrons of about 40 MeV for a 1 GeV primary proton, and 
the possibility of the reaction p + p — > n + + ... is about 1/3. 
Consequently, to maintain the annihilation flux at the observed 
level we should supply about 10 43 x 40 x ^ x 3 = 2.25 x 10 46 
MeV/s or 3.6 x 10 4l) erg/s. Since the time of positron thermal- 
ization is about 3 x 10 14 s we need the total energy in rela- 
tivistic protons ~ 10 55 erg, in order to produce the observed 
annihilation flux in the event of a single eruption. This amount 
of injection energy is possible by capturing a massive star with 
~ 30M G plus high conversion efficiency. However, the most 
natural explanation is a multiple capture scenario (multiple en- 
ergy injection model). Since the capture rate is ~ 10~ 5 yr , 
there would be about 100 captures in the past 10 7 years and 
each of them only requires ~ 10 53 erg, which is the most typ- 
ical value in Eq. (5). This non-stationary model also gives a 
natural explanation as to why the gamma-ray intensity emitted 
from the central high density region is so low. 

In Fig. 1 , we can see that both annihilation flux and gamma- 
ray flux in large regions are almost constant whereas the 
gamma-ray flux in the central high density region is much more 
sensitive to time. Furthermore the energy requirement of the 
multi-injection model is much less than that of the single cap- 
ture model due to accumulation of positrons in the thermal en- 
ergy region (Fig. 2). Our calculations show that this release 
should be about5xl0 52 erg i.e. a capture of one solar mass stars 
once a period t, ~ 10 5 years which is quite enough to generate 
the observed annihilation emission. This value is much smaller 
than the energy release in primary protons required for a single 
capture model (Cheng et al. 2006) which is about 10 55 erg for 
the intercloud gas density ~ 1 cm -3 . 




Fig. 1. The time variations of annihilation emission, gamma-ray flux from the central high density region (50 pc), and the flux of 
gamma-rays produced by protons escaping from the central core into the low density region (500 pc). 



4.2. MeV-GeV gamma-rays 

In Fig. 3, we calculated the wide band spectrum. In this pa- 
per we compa re the model calculations with t he observed data 
of EGRET dMaver-Hasselwander etalJI 1998b and COMPTEL 
dStrong et al.ll2005l) . In IStrong et all d2005b they summarized 
the data of EGRET from 30 MeV-30 GeV, COMPTEL from 
1 - 10 MeV and INTEGRAL from 100 KeV-0.5 MeV for the 
more extended region. The input energy in protons is chosen to 
reproduce the observed annihilation line flux. We can see that 
the model curve is consistent with the EGRET data in 5° radius 
but is substantially lower than the COMPTEL, which means 
that the in-flight annihilation of our multi-capture model cannot 
produce the observed 1-10 MeV photons whereas the neutral 
pions decay is capab le of explaining the GeV p h oton emission . 
Based on results of iBeacom & Yiiksell (12006b. iTotanil d2006h 
suggested that the model proposed by lCheng et al. I d2006b" will 
produce too much 1-10 MeV gamma-rays through in-flight 
annihilation and hence does not satisfy the observed data of 
COMPTEL. However, with the detail calculations by solving 
the proper kinetic equation, we can show that the in-flight 
annihilation will not over produce 1 -10 MeV. On the con- 
trary, calculations in the framework of lCheng et al.l (|2006) can- 
not produce enough 1-10 MeV gamma-rays to explain the 
COMPTEL data as shown in Fig. 3. We want to emphasize that 
most positrons are produced in the central high density region 
(~ 50 pc). Then they take ~ 10 7 years in the low density re- 
gion (~ 500 pc) to become thermalized positrons, which can 



capture electrons to form positronium. Our model injection of 
positrons in the low density region is equivalent to a constant 
injection of positrons with energy lower th an 10 MeV, which 
does n ot violate the constraint concluded by [Beacom & Yiiksell 
d2006h . 

However, in the energy range above 10 MeV, the in-flight 
annihilation emission is quite significant and gives the main 
contribution at ~ 30 MeV. It is interesting to note that at 
this energy the in-flight flux exceeds the contribution from 
inverse Compton scattering of relativistic electrons which is 
essential in the Galaxy especially at relative ly high latitudes 
(Dog iel & Ginzburgl 19891: (Strong et al.ll2000b . 



5. Discussion and conclusion 

We suggest that the observed diffuse GeV gamma-rays and the 
511 keV annihilation flux from 5° within the GC are conse- 
quences of multiple stellar captures by the Galactic black hole. 
The average injection energy carried away by relativistic pro- 
tons is ~ 6 x 10 52 erg per capture every 10 5 years. Such energy 
of the injection rate is sufficient to explain the observed elec- 
tron/positron annihilation flux. In our model, most positrons 
are produced in the central high density region (50 pc), and dif- 
fuse through the low density region to a distance ~ 500 pc dur- 
ing their thermalization time. Every capture takes place once 
every 10 5 years. Therefore the observed positron annihilation 
emission results from a population of thermalized positrons, 



Cheng et al.: Diffuse Gamma-ray Emission 



7 




10" 



10 u 



10' 
E (keV) 



10 



Fig. 2. The steady state distribution of positrons in the low density region (500 pc). 



which are produced, cooled down and accumulated by past 
hundreds of capture events instead of a single injection process. 
The gamma-ray intensity emitted from the central 0.5° region 
is very sensitive in time and the current intensity is most likely 
much weaker than its maximum value. Our numerical calcu- 
lations indicate that the in-flight annihilation cannot produce 
enough 1-10 MeV gamma-rays to explain the observed data by 
COMPTEL. On the other hand, unlike the diffuse gamma-rays 
in GeV range, the diffuse 1-10 MeV gamma-rays do not have 
a notable concentration within 500 pc. It is possible that they 
have different origins. 
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